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PREFACE. 



This paper presents the results of isentropic analysis in 
depicting oceanic flow patterns in a region of 'enow upwelling. 

Undertaken as a requirenent for the degree of has ter of Science 
in Aerology, this investigation was conducted at the U. S. haval 
Postgraduate School, Lonterey, California during the period fron 
August 1952 to January 1953. 

Grateful acknowledgment is made for the advice and assistance 
rendered by Assistant Professor Jacob D. V/ickhan in the preparation 
of this paper. 
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I. INTRODUCTION 



The phenomenon of upwelling, present wherever surface divergence 
exists, is of special interest in coastal regions. It is of profound 
commercial importance locally, as it pertains to the life cycle of the 
California sardine. Indeed the commercial importance of the sardine 
led to the organization of California Cooperative Sardine Research 
Program to investigate the virtual disappearance of the species from 
commercially accessible fishing grounds. The oceanographic observations 
collected in the routine operations of this program form the data basis 
for this paper. Aside from the commercial aspects of the upwelling 
process, it is of importance to naval operations in ocean areas. 
Successful submarine detection depends, to a large degree, upon the 
vertical temperature distribution in the sea which is profoundly 
effected by the process of upwelling. 

A qualitative explanation of the process based upon Ekman’s theory 
of wind driven currents was first suggested by Thorade 03 in 1509 and 
later developed by KcEwen Cl] in 1912. Elcnan’s theory, in part, states 
that the mass transport caused by wind stress upon the sea surface is 
directed at right angles cun sole to the wind direction. The wird acts 
to the depth of effective frictional resistance. Ekman [8, p 494] 

derived the empirical relation D = ^6 for this depth. It can 

fsvv* 

be seen from this equation that the depth to which the influence of the 
wind is felt depends upon latitude ar.d wind speed. The mass transport 
as given by T = K^ (8> P 49b] is a function of latitude and 
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t*. in turn is a function of V* as can be seen by the expression for 
stress ^8, p 490} where V is the wind at 15 meters 

above the surface. 

Consider a coast in the northern hemisphere with a wind blowing 
parallel to the coast such that the coast lies to the left of the 
wind. The mass transport of the wind driven current will be confined 
to the levels above O and will be directed away from the coast, thus 
the surface waters will be transported off shore at a rate proportional 
to the square of the wind speed. The continuity of the system requires 
that sub-surface waters rise to replace this mass deficiency. During 
the spring and summer months the semi -permanent Pacific anticyclone is 
established providing fairly steady winds roughly parallel to the 
central and southern California coasts, and it is during this period 
that upwelling becomes a conspicious feature. 

The phenomenon has been widely examined since the qualitative 
explanation was developed. Special aspects of the problem have re- 
ceived attention with emphasis upon the velocity of the vertical 
notion £ 2*] and the depth from which water is brought to the surface 1 7"] 
A more detailed investigation of the phenomenon was undertaken by 
Sverdrup ^ 6"] in 1938 using a single station line at right angles to the 
California coast off Port San Luis. From three scries of observations 
taken in 1 larch, Kay and July of 1937 along this station line, Sverdrup 
was able to develop a vertical section depicting the circulation 
perpendicular to the coast. The measured vertical motion of sub- 
stantial surfaces coupled with the principals of the equation of 
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of continuity to obtain horizontal components of otion was the 
analysis tool. This investigation, however, was based upon United 
data and assumed uniform velocity parallel to the coast. I cEwen ^3} 
has investigated the horizontal circulation associated with large 
scale horizontal eddys forming in conjunction with the upwelling 
process using the dynamic concept of the conservation of angular 
momentum. Fowever, no investigation of the regions of upwelling 
have, to the author’s knowledge, been made upon the basis of 
isentropic analysis. 4 

The validity of flow patterns inferred from isentropic analysis 
has been discussed by Parr [5] and Montgomery t4]. A general review 
of the concepts of isentropic analysis seems appropriate before pro- 
ceeding further. As is well known, surfaces in the atmosphere defined 
by constant potential temperature are surfaces along which adiabatic 
interchange of air masses, assuming air to be a perfect gas, result in 
no change in the distribution of mass. These surfaces are then isen- 
tropic surfaces. The question arises, are there surfaces in the sea 
which are isentropic. 

Montgomery ^4 J states that 

there exists no patent mechanism for altering the potential 
density of any water particle below the layer of direct surface 
influence. Therefore no flow of major proportions can take 
place across surfaces of constant potential density. 

Actually, as demonstrated by Sverdrup t®. P 417] , there exist no 

true isentropic surfaces in the sea; however, potential density 

surfaces are quasi-isentropic, implying that along these surfaces 

interchange of water masses result in only small changes in the mass 
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distribution. It is possible to deterr.'dne potential density surfaces 
in the sea. however, the normal density parameter observed in ocean- 
ographic work is , where (Ft is defined asCj^- \ooo . ^ 

the density at temperatures and salinities in situ but at atmospheric 
pressure. The difference between and , the potential density, 
is negligible above 1000 neters. This is demonstrated by Montgomery £ 4 J , 
who shows that the adiabatic temperature changes involved in bringing a 
water particle from depths less than 1000 neters to the surface arc small. 
Thus it is permissible to use Ct. surfaces as quasi-isentropic surfaces 
in water of depths less than 1000 neters. 

As in atmospheric isentropic analysis, the distribution of a con- 
servative property on the surface is employed in making inferences as 
to the flow pattern. Moisture concentrations , as shown by specific humidity, 
are used as the conservative property in the atmosphere; in the ocean 
salinity is the conservative property used. If a maximum of the con- 
servative property is found on the isentropic surface, it will, under 
the influence of a current on the surface, move with the current. As it 
moves, isentropic mixing along the surface will tend to reduce the in- 
tensity of the maximum in the direction of the current. If, at the sane 
tine, high concentrations are maintained at the original point of the 
maximum by some other process, then a tongue of maximum, but decreasing, 
concentration indicates a current along the tongue axis with flow 
toward low concentrations. Although there are no true isentropic sur- 
faces in the oceans, lateral nixing along surfaces is quasi-isentropic. 



Waters of differing salinities nixing along these surfaces will 
retain, approximately, their original densities; i.e., the products 
of nixing will remain upon the surface. 

Consider a concentration of high salinity water upon a 0^ 
surface, with a source of high salinity water to maintain the maximum. 
If a current on the surface moves tlrrough the area of maximum salinity, 
a tongue of high salinity will extend in the direction of current flow 
in a manner similar to the situation shown in atmospheric isentropic 
analysis. Thus flow patterns can be inferred from the distribution of 
high or low salinity tongues on surfaces. This is, of course, a 
generalization of Wust’s I'ernschicht method of tracing absolute 
salinity naxir.uns in the vertical, as show; by hontgonery ^ 4 ~\ . 

It is apparent from this discussion that isentropic analysis as 
applied to the oceans is primarily qualitative rather than qualitative. 
The descriptive nature of such analysis is, of course, an end in itself 
and can provide valuable info mat ion without the knowledge of actual 
current velocities. Possible quantitative results will be discussed 
in a later chapter. 

The use of salinity as a conservative property nay be questioned 
on the grounds that density is not independent of salinity; however, 
density does not depend entirely upon salinity. If there existed a 
unique relationship between density and salinity then salinity sur- 
faces and density surfaces in the ocean would be parallel. However, 
these surfaces do intersect, so it would seen that the distribution of 
salinity upon a 0^ surface is brought about by the flow patterns and 
the raxing which exist. 
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II. DATA USED IE' INVESTIGATION 



The data used in this investigation were provided by the operations 
of the California Cooperative Sardine Research Program. The raw data 
has been processed and published by Scripps Institute of Oceanography. 

The area investigated is shown in Plate I which includes the positions 
and number designators of the oceanographic station used. These stations, 
in the course of the program, were occupied monthly by vessels of the 
Research Program and observations of the following parameters were made: 

(1) T 

(2) S 

(3) O* 

(4) P0 4 -P 

The values of these were interpolated graphically to standard depths, and, 
from (1) and (2), values of the following parameters were determined for 

standard depths: 

(5) <3\ 

( 6 ) \O s <S 

(7) 

The parameters used in this investigation were: 

(1) D 

( 2 ) S 

It was necessary to assume that the data were synoptic in nature. 
Actually the tine interval involved in completing the observations for 
the area under investigation varied from a minimum of 5 days in Narch, 



( 6 ) 



when not all stations were occupied, to a maxima* of thirteen days 
in Kay. Since conditions in the oceans change slowly, this assumption 
does not invalidate the resulting general features of the flow; but it 
does make any interpretation of the niror features doubtful. 

It is also necessary to assume that station positions were fixed 
from month to month. This assumption is not quite fulfilled, but the 
error is snail since the maximum deviation from one month to the next 
was two nautical miles. The error involved in this assumption is 
probably negligible compared to errors in the assumption of synopticity. 

All parameters were assumed to vary linearly between standard depths. 
This assumption reduced the labor and time of interpolation considerably. 
Since the data as received was interpolated to standard depths, it was 
felt that little accuracy could be gained by vertical plots of the para- 
meters in an attempt to reproduce the original sounding. In view of the 
above, and because the results are intended to be qualitative, the tine 
saving computational method of interpolation was felt to be justified. 

The file of analyzed weather charts of North America and Pacific 
ocean areas maintained by the Naval Postgraduate School, Konterey, 
California, were utilized to determine wind components parallel to the 
coast during the period investigated. The calculated values of the 
wind component are qualitative in nature. This is true because the 
area concerned represents a small segment of the analyzed weather chart 
and, in many cases, isobaric spacing could be adjusted through wide 
ranges, while still satisfying the surface reports. In all cases 
possible, results of geostrophic measurements were checked against 
observed winds to reduce errors. 
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III. IIETiiOD CF AhALYSIS 



As a first step, the conponent parallel to the coast of the 
average wind over the area was computed in this ratter. The average 
geos trophic wind over the area was measured directly from, the daily 
1230 GOT analyzed chart for the period from 1 February to 31 July. 

The winds were assumed to deviate ten degrees from the isobars, and 
corrections for isobaric curvature and air rass stability were made 
in accordance with the procedure outlined in II. 0. 604 [_10, pp 16-17} . 

An average stability factor of .60 was selected for the entire period. 

This represents a sea temperature minus air temperature difference of 
from one to seven degrees. The components of these winds parallel to 
the coast were computed and five day means of these values were plotted. 
The results are shown in Plate II. The hatched areas in the figure 
represent the time interval involved in the cruises during each month. 
Since the stress exerted by the wind upon the sea surface is a function 
of the square of the wind speed, it would seem appropriate to graph these 
values rather than the actual wind speeds, but in this investigation only 
the qualitative aspects of the wind in producing upwelling is necessary. 

Inspection of the diagram shows considerable variability in the 
coastal components, but several features are immediately apparent. The 
winds conducive to upwelling begin early in February and continue through- 
out the period. The period between hay and June show the weakest and 
most erratic values, while the eighteen day periods prior to the hay and 
July oceanographic observations shew relatively steady components of 
moderate velocities following periods of maximum velocities. From this 
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it is inferred. that tlie beginnings cf the upwelli ig si oulcl be appare: t 
in I arch ard that the r.axi! um upwellir.g should be apparent in I ay and 

My. 

The second step i' the investigation was the preparatio; cf 
isentropic charts for the months of 'ebruary through July. he 
surf?vCes 0"t~ 2®*00 ai 10^*26. SO were selected except for the month 
July, in which south the additional surfaces O'*. -26.75 2 . 1-00 

and 27-iS were also investigated. his additional analysis was 
uudertahen in order to obtain a more detailed picture in July which, 
on the basis of thc<5t - 26 .oo so 

be the north of most pronounced upweliii g (see Chapter IV). 'he dis- 
tribution of sali ity on these selected surfaces was deter' -ired by 
linear interpolation fro;.: values at standard depths. Also the geometric 
depths of the selected surfaces were siilarly deterr.ined. hese two 
values, salinity on the surface and the depth of the surface, were 
determined for each station for each month. Isohalines were then 
drawn on each surface; and, on the same chart, the bathyr. etric topography 
of the given surface was also dram. The results are shorn, in Hates 

iii-xvn . 
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i”. i: rmp;r.:?A;ic. of uiiULfb a..i cu.olusiclj 

As discussed in Chapter I, salir.it;' distributions c.. <Tt surfaces 
are assumed to indicate flow patterrs. In accordance with this, arrows 
depicting the general flow pattern have beer, constructed on the basic 
charts (Plates III-XVTI). Inspection of the charts shows two pror .inert 
features, first, a large cyclonic eddy is evident or. all charts except 
the 27.25 surface for „uly, where the pattern breaks doi/n at these 

great depths. This eddy’s presence is inferred from the tongue of low 
salinity which extends toward the coast in the southern portion of the 
area and then curves northward inshore. Second, the isolated naximm s 
of salinity occur on don.es in theOt surfaces, wliile conversely isolated 
r.dnir.ums occur on depressions in the surfaces. The first feature is not 
surprising since the shoreward extension of low salinity water and the 
existence of the large cyclonic eddy have been demonstrated by other 
methods in previous investigations of the sane general area £ll, p 2?'] . 
This feature then, as depicted by iser.tropic aralysis, substantiates 
other findings. The second feature is of wore interest to the author. 
Under the assumption thatO^ surfaces are substantial surfaces it 
follows that, if vertical notion occurs within a region, the surfaces 
Must move with the water. It will be recalled that, on the basis of 
the longshore wind component calculations (see Plate II), Larch was the 
month during which upwelling should first occur, figure 1 demonstrates 
the effect of the upwelling upon theQ^ surfaces during the irterval 
between February and Larch. The vertical section shown in Figure 1 is 
for station line 80. Upward notion is indicated inshore of station 80-70 
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Vertical section through Station Line 80 showing 
vertical displacement of selected <1* surfaces 
fro® February to f'arch 1950 

February surfaces 

f'arch surfaces 

depth in reters 

Figure 1 
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to depths of at least 500 meters, while downward motion is shown 
west of station 80-70 to depths below 200 r.eters. This downward motion 
is to be expected west of the area of strongest longshore wind components 
since the decreasing wind stress would produce decreasing offshore trans- 
port and surface convergence would occur, necessitating notion do’./nward . 
/he decreasing wind velocity and, thus, wind stress is to be expected 
since the pressure gradient in the high pressure cell giving rise to the 
winds decreases toward the center of the high. 

The average upward displacement of the CJt surfaces in the inshore 
area is approximately 40 i eters, and a sir pie calculation gives the 
average vertical velocity, based upon the surface displace :e:~ts, as 

1.5 x 10“3 cr’/sec. If the surface divergerce rcspoi sible for this 
upwelling continues, notion must tahe place through the CTt surfaces 
because these surfaces seen to become stationary near, but below, the 
sea surface. This Leans that, in the region of continuing surface 
divergence, the compensating inflow from below cannot be isentropic, 
since, by definition, any component of motion which crosses the <3^ 
surfaces is non-isentropic. If the cor ponent of the wind parallel to 
the coast can be used as an indicator of surface divergence, then 
surface divergence was a continuous phenomenon tluoughout the period 
under investigation;. Inspection of the data for this period shows that 
salinity increases with depth at any given station except for r iuor 
variations. This means that no major salinity maximums or mini: .urns 
exist in the vertical. This observation offers a clue to the origin 
of the salinity r aximums and ninimums appearing on the CEt surfaces. 
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namely, that water characterized by high salinity on a given 
surface cones from depths below the level of that surface. Conversely, 
low salinity water on a given surface cores fror depths above the level 
of that surface. If the above reasoning is valid, then salinity ex- 
tremes on a given surface indicate vertical i ixing, which is non- 
isentropic, or some other ron-isentropic process. 

hontgomery^ 4} has discussed a non-isentropic process to allow for 
the upwelling necessary to compensate for surface divergence. This 
process he calls nixing upward and concludes that it is due to vertical 
turbulence occurring at the boundary of the homogenous layer. He further 
concludes that this vertical mixing cannot extend deeper than the lower 
boundary of the layer of surface influence (see Chapter I). Assuring 
that mixing upward will compensate for the surface divergence, it becomes 
necessary to explain the fonntion of the salinity maximum on the elevated 
Ct stirfaces near the lower boundary of the layer at surface influence. 

An example of this situation is shown on Plate XIII where theQTt 5 26.00 
surface lies between 50 and 75 meters in the region southwest of the 
Channel Islands with the associated salinity maximum on the C3t. surface. 

A vertical temperature sounding for Station 83-60 which lies in the region 
of the elevated surface, is shown in Figure 2 with the depth of the 
- 26.00 surface indicated on the sounding. The figure shows that the 

OTt - 26.00 surface lies at the lower boundary of the layer of surface 
influence. Consider a point on this boundary; if water from below is to 
penetrate to levels above this boundary, then the density of this pene- 
trating water nust be altered by some non-isentropic process. If this 
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were rot true tLe OE = 26.00 surface would rove upward ard penetrate 
to the surface , and this is rot observed. According to the observation, 
the local rate of change of density, ~~ , at tl e lower boundary of 

O t. 

the surface layer rust equal 0. he rate of char 0 e of density of a 
roving particle nay be expressed as the sun of the local and advective 

rates of change of density: 

- iLe + a + ^ s^ur 

dkt <=>•* 

Assuming horizontal advection to be small in the region of maximum 
upward mixing, this expression reduces to 

= ur 

The individual rate of change of density is, under these assumptions, equal 
to the vertical advection of density. A possible explanation for the de- 
crease in density of the water particle as it moves into the surface layer 
is by eddy diffusion of heat from the warmer surface layer. This eddy 
diffusion, or mixing, does not appreciably decrease the salinity, however, 
as salinity gradients in the nixed layer are very small. The heat content 
of this layer is replenished, in turn., by solar radiation. Thus the 
density of the particle could be decreased without appreciably decreasing 
the salinity, resulting in a concentration of high salinity water appear- 
ing on the elevated Q\ surface iixediately below the surface layer. 

It is interesting to note that calculations, >‘ot preserted here, using 
the concept of density decrease by solar radiation ^ive a vertical 
velocity, of the upward moving particle, of the sar e order 






( 15 ) 



of nagnitude as the vertical velocities of upward notion deduced from 
the vertical displacement of the <St surfaces early in the upwelling 
period. 

The salinity maxiruins evident upon thefts 26.00 surfaces, when 
they are close to the sea surface in other months, can be explained by 
the same reasoning, if this reasoning is valid. 

The existence of salinity minir.uns on the deeper surfaces and 
the existence of both maximums and mininums at intermediate depths, 
which are particularly evident on theQ^'a 26.75 and r 27.00 
surfaces for July (see Plates XV and XVI), are probably due to vertical 
mixing, but the cause of such mixing is not readily apparent. Internal 
waves, characterized by large velocity shear across interfaces, are a 
possible general explanation of pronounced vertical mixing at depths 
where salinity maximums and riininuns appear. Another factor to be 
considered is the effect of bottom friction in producing turbulent 
vertical mixing. However, it is beyond the scope of this investigation 
to assess the influence of these factors upon the observed salinity 
distributions . 

As discussed in Chapter I, the concept of upwelling associated 
with wind driven currents requires a transport of water perpendicular 
to the wind direction. In the area under investigation, this transport 
is, of course, directed offshore, evidence of this offshore transport 
is inferred from the salinity distribution on the ~ 26.00 surface 
for Kay (Plate IX) and the « 26.00 surface for July (Plate XIII). 
The pronounced high salinity tongue extending offshore from the salinity 
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maxinur.s in the region south of the Channel Islands is, in each case, 
considered to give evidence of this offshore transport. As discussed 
in Chapter III, the results of the longshore wind component computations 
indicated that hay and July should he the months of maxi] un upwelling, 
and the observed salinity distribution on the upper surface indicates 
that the offshore transport was -ost pronounced during these months. 

Longshore transport of the more saline upwelled water is inferred 
from the southward extension of high salinity tongues on the = 26.00, 
26.50, 26.75 and 27.00 surfaces for July (Plates XIII-XVI). This result 
is in qualitative agreement with previous investigations of the area 
111, P 27} which indicate that the more saline waters are carried south 
from the region of upwelling by the eastern edge of the California 
current . 

The results of this investigation have verified, in the author 1 s 
opinion, results previously obtained by different methods; and, as such, 
have proved of some value. Lore investigation of the processes giving 
rise to the observed salinity maxixurs and minirums or surfaces is 
necessary before these distributions can be said to be of significance 
in picturing the centers of upwelling or subsidence. The qualitative 
nature of isentropic analysis can be overcome to some extent by con- 
siderations of diffusion and advection as effecting the distributions 
of so r e conservative property, e.g., salinity. Such calculations were 
beyond the scope of this investigation. 
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Plot of wind component parallel to coast, in knots. 

1 February to 31 July 1951. Shaded areas represent 
tine involved in completing oceanographic observations. 
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